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yA  tractable,  five-parameter  family  of  continuous,  unlmodal  probability 
distributions  Is  developed.  Special  cases  Include  the  Bernoulli  trial,  uniform,  power 
series,  exponential,  triangular  and  Laplace  (double  exponential)  distributions.  The 
family  has  closed-form  density  function,  cumulative  distribution  function.  Inverse 
cumulative  distribution  function,  hazard  function,  and  residual  moments.  Statistical 
properties,  parameter  determination  and  random  variate  generation  are  discussed.  Two 
numerical  examples  are  given. 


KEY  WORDS:  Input  modeling}  Monte  Carlo^ Process  generation,  Simulation,  Statistical 


modeling. 


1.  INTRODUCTION 


We  consider  the  five-parameter  family  of  probability  distributions  having  density 
function 


[1  +  (x-a)(6/e)1/2/(01)]^^s  /  (/W1/2) 


f(x)  = 


[l  -  (x-a)(e/8  )1/2/(^(l^7»](w)/f  /  (0(6e)l/2) 


0 


If  a  <  x  <  a 
If  a  <  x  <  b  (l) 
otherwise 


where  a  =  a  -  0  7  ( t/8  )l/ 2  Is  the  lower  bound  and  b  =  a  +  0  (1-7)  (6/e),/2  Is  the 
upper  bound  of  the  random  variable  X.  The  parameters  have  ranges  -  00  <  a  <  00, 
0  <  0  <  00,  0<7<1,0<6<oo,  and  0  <  e  <  00.  This  family  can  model  a  wide 
variety  of  univariate,  continuous,  unlmodal  probabilistic  phenomena  by  varying  the  five 
parameters.  Properties  of  the  family  are  developed  In  Section  2,  parameter 
determination  Is  discussed  In  Section  3,  and  two  numerical  examples  are  given  In 
Section  4.  The  rest  of  this  section  Is  a  brief  review  of  related  literature. 

Schmelser  (1977)  surveys  the  literature  of  versatile  families  of  distributions. 
Including  the  Pearson  and  Johnson  systems  of  distributions;  approximations  to  the 
Inverse  distribution  function  via  expansion  techniques,  polynomlnal  regression,  and 
rectangular  approximation;  and  various  four  parameter  distributions.  Johnson  and 
Kotz  (1970)  review  various  classical  methods  In  some  detail.  We  focus  here  on  the 
more  recent  literature,  for  which  the  motivation  often  has  been  to  develop  families 
useful  as  Input  models  In  Monte  Carlo  simulation. 

Johnson  et  al.  (1980)  develop  a  four-parameter  family  of  symmetric  distributions 
that  Includes  the  uniform,  normal,  and  Laplace  distributions  as  special  cases.  Since 
random  variate  generation  Is  a  simple  transformation  of  a  gamma  random  variate  and 
since  the  kurtosls  Is  a  simple  function  of  the  parameters,  the  family  Is  Ideal  for  Monte 
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Carlo  studies  of  the  effects  of  departures  from  normality.  Fitting  the  distribution  to 
data  Is  not  discussed. 

Stacy  (1962),  Stacy  and  Mihram  (1965),  and  Harter  (1967)  discuss  the  four- 
parameter  gamma  distribution,  which  Is  based  on  taking  a  gamma  random  variable  to  a 
power.  Both  the  three- parameter  Welbull  and  gamma  distributions  are  special  cases. 
This  distribution  Is  useful  for  modeling  lifetimes,  although  parameter  estimation  Is  not 
straightforward. 

Hora  (1983)  develops  a  modeling  procedure  based  on  the  Inverse  distribution 
function.  The  analyst  selects  a  tractable  distribution  function  F  that  Is  a  reasonable 
approximation.  The  model  Is  F-1(g(p)),  where  g(p)  Is  a  polynomial  model  fitted  by 
regression  to  the  data.  Random  variate  generation  and  data  fitting  are  straightforward, 
but  the  analyst  must  provide  the  distribution  function  F.  The  families  discussed  below 
provide  good  candidates  for  F. 

Burr  (1942,  1973)  developed  a  more  flexible  distribution  with  the  same  range  and 
straightforward  random  variate  generation.  Heavy  and  light  tailed  distributions  are  not 
obtainable,  nor  are  symmetric  distributions,  but  the  family  Is  useful  for  modeling 
lifetimes,  although  again  parameter  estimation  Is  difficult. 

Ramberg  and  Schmelser  (1972,  1974)  generalized  Tukey’s  (1960)  lambda 

distribution  to  obtain  a  family  for  all  but  light  tailed  distributions.  Including  symmetric 
distributions.  The  exponential  distribution  Is  a  limiting  case  and  the  normal 
distribution  function  Is  approximated  within  a  tolerance  of  .001  In  the  distribution 
function.  (See  Joiner  and  Rosenblatt  (1971)  for  a  related  approximation.)  Random 
variate  generation  Is  straightforward.  Dudewlcz  et  al.  (1974)  and  Ramberg  et  al.  (1979) 
provide  tables  to  fit  the  distribution  by  matching  moments.  Other  estimation 
techniques,  such  as  maximum  likelihood  estimation,  are  not  straightforward  due  to  the 


distribution  function  and  density  function  not  being  expressible  In  closed  form, 
•• 

although  Ozturk  and  Dale  (1985)  consider  minimizing  the  squared  errors  between  the 
empirical  and  fitted  distribution  functions. 

Schmelser  and  Deutsch  (1977)  describe  a  family  of  distributions  that  can  be  used 
to  obtain  a  distribution  having  any  first  four  moments.  The  exponential  distribution  Is 
a  limiting  case.  Closed  form  expressions  are  given  to  calculate  parameter  values  to 
match  the  desired  mode,  a  percentile  of  the  mode  and  any  other  two  quantiles. 
Moments  may  be  matched  graphically  or  via  nonlinear  programming.  Random  variate 
generation  Is  closed  form  and  requires  only  one  exponentiation  operation.  Kottas  and 
Lau  (1979)  use  the  distribution  In  Inventory  modeling  as  the  distribution  of  lead-time 
demand  and  Lau  and  Zakl  (1982)  use  the  distribution  to  test  sensitivity  to  lead-time- 
demand  assumptions.  The  disadvantage  of  this  distribution  Is  that  while  any  four 
Independent  properties  may  be  specified,  the  shape  of  the  distribution  Is  only 
satisfactory  as  a  rough  approximation.  The  shortcomings  Involve  truncated  tails  and  a 
density  function  value  at  the  mode  that  Is  often  either  zero  or  Infinity.  This  family 
could  be  generalized,  with  little  loss  of  tractablllty,  to  five  parameters  analogous  to  the 
family  developed  here. 

The  family  of  distributions  developed  here  has  continuous  densities,  obtains  any 
mean  and  variance  as  well  as  a  wide  range  of  commonly  used  third  and  fourth 
moments,  has  straightforward  random  variate  generation,  has  closed-form  moments, 
and  has  tractable  parameter  determination. 
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2.  PROPERTIES  OF  THE  FAMILY 


As  follows  directly  from  the  density  function,  the  distribution  function  Is 


0 


F(x)  = 


^[l  +  (x-a)(6/e)1/V(/?'Y)]1/‘ 

1  -  {(l-7)[l  -  ( x-a)(e/8  f'2  /  (£  (1— ^y))]1/c> 


If  -  oo  <  x  <  a 
If  a  <  x  <a 
If  a  <  x  <  b 
If  b  <  x  <  oo 


where  again  a  =  a  -  [3  7  (e/8)l/2  and  b  *=  a  +  /?  (1  -  7)  (6/01/2  are  the  lower  and  upper 
bounds. 

An  appealing  feature  of  this  family  of  distributions  Is  that  the  role  of  each  of  the 
five  parameters  Is  relatively  Intuitive  and  Independent  of  the  others.  The  location 
parameter,  a,  Is  the  mode  or  antimode.  The  scale  parameter,  /?,  Is  the  range  when 
8  —  e  and  Is  always  proportional  to  the  standard  deviation,  which  Is  also  a  function  of 
7,  8,  and  €.  These  three  remaining  parameters  determine  shape,  but  again  each  has  a 
specific  role.  The  probability  to  the  left  of  the  mode,  a  rough  measure  of  skewness,  Is 
7;  l.e.,  P{X  <a}  =  7.  The  shape  of  the  distribution  Is  determined  primarily  by  8  to 
the  left  of  a  and  by  c  to  the  right  of  a.  The  shape  Is  symmetric  If  and  only  If  5  =  e 
and  7  =  1/2.  For  8  >  1  (e  >  1),  the  curve  to  the  left  (right)  of  the  mode  will  be  J- 
shaped  with  skew  to  the  right  (left).  For  6  <  1  (e  <  1),  the  curve  to  the  left  (right)  of 
the  mode  will  be  J-shaped  with  skew  to  the  left  (right). 

These  parameters  yield  a  wide  range  of  shapes,  as  can  be  seen  by  considering  some 
special  cases.  First  consider  symmetric  distributions,  for  which  8  =  e  and  7  =  1/2.  As 
8  —  e  approaches  Infinity  with  o  —  7  =  1/2  and  0  =  1,  P{X  =  0}  — ►  1/2  and 
p{X  =  1}  — ►  1/2,  a  fair  Bernoulli  trial.  The  uniform  distribution  over  the  Interval 
(a,b]  Is  obtained  by  setting  a  =  a  +  (b-a)  7,  0  =  b-a,  8  =  e  =  1.  (In  this  one  case,  7 
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may  assume  any  value  In  the  unit  Interval  and  the  distribution  remains  symmetric.)  As 
6  =  e  approaches  zero  with  7  —  1/2,  the  limiting  distribution  Is  the  Laplace.  Thus  for 
symmetric  distributions  the  family  provides  a  kurtosls  range  from  one  (Bernoulli  trial) 
to  six  (Laplace).  Very  heavy  tailed  distributions  cannot  be  obtained. 

A  wide  range  of  asymmetric  shapes  Is  also  possible.  As  6  —  e  approaches  Infinity 
with  a  —  7  and  0—1,  P{X  =  0}  — ►  7  and  P{X  =  1}  — *•  1-7,  a  biased  Bernoulli  trial. 
By  letting  a  —  0,  0  =  X/e,  7  =  0,  and  6  =  e,  the  limiting  case  as  e  -*■  0  Is  exponential 
with  mean  X,  as  can  be  seen  from  by  considering  nonnegative  values  of  x  In 

11m  F(x)  =  11m  1  -  [1  -  x/(X/e)]1/£  =  1  -  e_x/x 

(  ->  0  f  — »  0 

which  Is  the  distribution  function  of  the  exponential  distribution  having  mean  X.  The 
power  series  distribution  having  density  function 

X  xX-1  for  0  <  x  <  1 

f(x)  = 

0  elsewhere 

Is  obtained  when  a  —  0  —  7  =  1  and  8  =  e  =  1/X. 

Some  graphs  Illustrate  the  versatility  of  the  family.  Figure  1  shows  graphs  of  the 
distribution  for  7  =  0.5  for  various  values  of  6  =  e.  Changing  a  and/or  0  would 
change  the  location  and/or  scale,  but  not  the  shape,  of  the  distribution.  Since  7  =  0.5, 
the  distribution  Is  symmetric  for  all  <5  =  e.  Figure  2  shows  graphs  corresponding  to 
'y  =  o.75.  Here  75%  of  the  area  under  each  curve  lies  to  the  left  of  the  mode  (anti¬ 
mode)  a  for  all  6  =  e.  When  6  =/=  e,  the  left  and  right  sides  of  the  distribution  differ, 
resulting  In  arbitrary  combinations  of  tall  shapes. 


Solving  for  x  as  a  function  of  p  In  the  distribution  function  yields  the  Inverse 
distribution  function 

a-  0  t(c/6  )1/2(l  -  (p/nr)tf)  If  0  <  p  <  Tf 

x  =  F'^P)  =  (2) 

a  +  0{l-^){6/ c)1/2[l-((l-p)/(l-Tf))<]  If  7  <  P  <  1 

which  gives  the  pth  quantile  of  X  as  a  function  of  p.  In  addition  to  calculating 
quantiles,  the  Inverse  transformation  Is  Important  as  a  method  of  random  variate 
generation.  If  P  ~  U(0,1),  then  F_1(P)  ~  X.  Random  variate  generation  based  on  (2) 
Is  relatively  fast  In  that  only  one  exponentiation  Is  required.  The  closed-form  Inverse 
distribution  function  allows  the  direct  generation  of  order  statistics  via  the  methods  of 
Lurie  and  Hartley  (1972),  Schmelser  (1978a),  Ramberg  and  Tadlkamalla  (1978),  and 
Schucany  (1972),  as  surveyed  In  Schmelser  (1078b).  The  Inverse  transformation  Is  also 
Important  In  simulation  experiments  as  a  method  to  Induce  positive  or  negative 
correlation  for  variance  reduction,  as  discussed  e.g.  In  Wilson  (1985). 

For  c  satisfying  F(c)  <  1,  the  residual  moments,  E{Xk  |  X  >  c},  when  a  =  0,  are 

w  k  (-1)'  (?)  (V*+1-p/{+1) 

E{Xk  |  X  >  c,  a  =  0}  =  p  j(-y)k  {e/6  £ -  ~6 - 


+  (l-y)k  {6/ef/2 


^  (-D‘  (*)  (i-puy<+i) 
So  (1«+1)  (1-7)'* 


/  [1  -  Pc) 


for  k  =  1,  2,  •  ■  *  ,  where  pc  =  F(c  ;  a  =  0),  pt  =  mln{7,  pc},  and  pu  =  max{7,  pc}. 
The  residual  moments  are  derived  In  the  Appendix. 


The  moments  centered  at  a  are  obtained  In  the  limit  as  c  — ►  -  oo,  In  which  case 
p*  =  0  and  pu  =  7,  yielding 

E{Xk  |  a  =  0}  =  f  {-  (-  '7)k+1  ( e/6  )k/2  [  E  (-l)‘  ( f )  /  0*  +  1)1 

i=o  s 

+  (l^r)k+1  «/£)k/2  (f)  /  (le+1)]} 

1=0 

=  p  {-■(-  l)k+1  (e/6  )k/2  B(k+l,  S'1) 

+  (l-^)k+1  (5/c)k/2  e~l  B(k+1,  e-1)} 

where  B(. , .)  Is  the  beta  function,  which  can  be  evaluated  numerically  or,  since  one 
argument  Is  Integer,  algebraically  using 

6-1  B(k+1,  6-1)  =  (*)  /  (\6  +  1) 

1=0 

==  p  ki  /  no$+i) 

i=o 

The  specific  four-parameter  case  of  6  =  e  yields 

E{Xk  |  a  =  0,  6  =  e}  =  p  ((i-'y)k+1  -  Hy)k+1]  S'1  B(k+l,  fi"1) 

In  the  general  five-parameter  case,  the  mean  Is 

E{X}  =  a  +  /?  (6e)1/2  [((l-/)2/(e+l))  -  tf/(6+l))) 


and  when  8  =  e 


E{X  |  6  =  e}  =  a  +  0  [(1-2 7)6  /  (5+1)) 


The  mean  Is  equal  to  the  mode  a  when  7  =  1/2  and  6  =  e  or  when 
7  =  [(6+1)  -  ((6+l)(e+l))1/2]/(6-e)  and  6  =/=  e.  Equivalently  the  mean  and  mode  are 
equal  If  t  —  [(1-7)2(6+1)/t2J  -  1  or  6  =  [t2(c+1)/(1-7)2]  -  1  . 

The  variance  can  be  calculated  using  V{X}  =  E{X2  |  a  =  0}  -  E2{X  |  a  —  0}. 
In  the  specific  four-parameter  case  of  6  —  e,  V{X  1 6  —  e}  = 
W  Sf  [27  (1-7)  (^-1)  +  1]  /  1(26+1)  (6+1)2]. 

3.  PARAMETER  DETERMINATION 

To  be  useful  In  statistical  modeling,  a  family  of  distributions  requires  methods  for 
determining  appropriate  parameter  values.  Two  contexts  —  (1)  obtaining  a  set  of 
specified  properties  such  as  moments  and  fractlles  and  (2)  fitting  to  data  —  are 
discussed  here.  The  mixed  case  of  obtaining  some  properties  and  fitting  to  data  Is 
discussed  briefly  In  Section  4. 

3.1  Obtaining  Specified  Properties 

Consider  fitting  a  distribution  to  a  set  of  specified  properties,  rather  than  to  data. 
Situations  In  which  no  data  are  available  often  occur  when  modelling  proposed  systems 
that  do  not  yet  exist.  The  set  of  properties  to  be  obtained  arises  from  technological 
considerations  and  from  expert  (or  nonexpert)  opinion.  The  former  Is  often  lower  and 
upper  bounds  and  the  latter  Is  often  low-order  moments. 

Four  sets  of  properties  are  discussed  here.  In  each  case,  appropriate  parameter 
values  can  be  found  using  closed-form  formulas,  sometimes  combined  with  a  single 
unldlmenslonal  search,  to  obtain  any  of  several  sets  of  properties.  For  each  set  of 


properties,  the  corresponding  parameter  values  are  given  below.  Their  derivation  Is 
straightforward  and  other  sets  can  be  derived  as  needed. 


Here  /i  and  o2  denote  the  mean  and  variance,  a  and  b  the  lower  and  upper  bounds, 
and  m  the  mode.  Also  pm  Is  F(m),  (plf  xt)  Implies  that  F(Xj)  =  pt  where  xx  <  a  and 
pj  <  7,  and  (p2,  x2)  Implies  that  F(x2)  =  p2  where  x2  >  a  and  p2  >  7. 

I.  Specify  a,  m,  pm,  (plf  xj,  and  (p2,  x2).  Then 

a  —  m 
1  =  Pm 

6  =  ln((xra)  /  (o:-a))/ln(p1/7) 

Solve  for  1  recursively,  by  trial  and  error,  or  any  unldlmenslonal  search 
procedure  In  e  =  8  [l  -  ((l-P2)/(l-nf))e]  ((I-^VtO  (<*-a)  /  (x2-a) 

0  =  ((a-a)/7)  (V01/2 

II.  Specify  a,  b,  m,  pin,  and  (p,.  Xj).  Then 

a  —  m 
•7  =  Pm 

6  =  ln((x1-a)/(a-a))  /  ln(Pi/7) 

€  =  <5  (a-a)  (1-7)  /  ((b-aft) 

0  =  ((a-a)/7)  (V01/2 

III.  Specify  a,  b,  /1,  and  m,  pm.  Then 

a  —  rn 
”7  ~  Pm 

0  =  [(a-a)  (b-  a)  /  (7  (1- -7))] 1/2 
c  —  0  7  /  (a-  a) 

Solve  for  e  In  the  quadratic  equation 

c2  [c(/i-a)  /  /?  +  72  --  c2(I-7)2]  +  e  [(/-*-<*)  (c2+l)  /  (/?c)  +  27-1} 


+  (fi-a)  /  (/5c)  =  0 


Then 


6  =  e  c2 

Here  c  =  (6/e)1/2  Is  an  Intermediate  value. 


Specify  a,  /+  a 2,  m,  pm.  Then 


a  =  m 


T  —  Pm 

Solve  for  8  using 

^  +  (#*  -  a)2  =  W(e/6)l(l/(26+l))  -  (2/(6+l))  +  1] 

+  (l  -  7)3(V£)[1  -  (2/(e+l))  +  (l/(2e+l))]} 

where 

e  =  — - 1 

6(a-a)  6+1 


0  =  ((a-a)/7)  (6/e)*/2 


The  four  sets  of  properties  are  certainly  not  adequate  for  every  situation.  But 
given  another  set  of  properties,  the  user  at  most  needs  to  solve  five  nonlinear  equations 
for  the  five  unknown  values  a,  0,  7,  8  and  e.  For  most  properties,  however,  algebraic 
simplification  reduces  the  problem  to  one  or  two  equations. 


3.2  Parameter  Estimation  from  Data 

The  more  classic  distribution  fitting  situation  Is  to  estimate  parameter  values  from 
data  Xj,  x2,  •  •  •  ,  xn  by  assuming  the  data  arise  Independently  from  some  particular 
member  of  the  family  of  distributions.  Common  methods  Include  maximizing 
likelihood,  matching  moments,  and  minimizing  a  goodness-of-flt  statistic  such  as  the 
Kolmogorov-Smlrnoff  statistic,  the  chi-square,  or  a  weighted  squared  deviation  between 
the  empirical  and  fitted  distribution  function.  Ad  hoc  methods,  as  discussed  In  the 
next  section,  for  fitting  the  five-parameter  family  work  reasonably  well,  but  a  well- 
developed  theory-based  methodology  has  not  yet  been  developed.  As  a  first  step,  the 
five  maximum  likelihood  equations,  obtained  by  setting  the  derivatives  of  the  likelihood 
function  function  to  zero,  are  stated  here. 

a  :  (1  -6  )Sh+  z  jtf/e)178]-1  =  (l-«)  £  [(1-7)  -  z,  (e/6  (3) 

i=\  i=n,+l 

n  (6  e)l/2  =  (6-1)  £  hr  +  h (V01/2]_1  Zj 
i=l 

-(€-1)  £  1(1-7)  -  Zj(e/6  )1/2]_1  Zj  (4) 

i=Di+l 

(1-0  £  Zj/((l-7)  -  Zj(e/5  )1/2) 

i=n!+l 

(1-6 )  Y,  Zi/(nr  +  Zj(6/e)I/2) 


7  :  (7-1  )/7  = 


(5) 


6:  6  =  (/?/n)  {  -  2  £  ln[l  +  (zt(6/0 1/2/7)] 

i=l 


+  d-s)  E  ft  +  *i(«/<)1/'T‘  2,(</0,/J 

i=l 

+  (1-0  E  1(1-1)  -,,(«/« )1/a]-I*i«/<),/s}  (0) 

i=n,+l 

c  :  e  =  (£/n)  {-  2  2  ‘“l1  “  zi(£/5  )1/2/(l-0r)l 

l==n,-n 

+  («-D  E  b  +  2i(«/€),/V  *i(e/«  ),/2 

i=l 


+  (*-i)  e  [(i-r)-2,(</4),/!r'il(</'5  ),/a>  p) 

i*=  nj+l 

where  Zj  =  (x,-a  )/0  and  nt  Is  the  number  of  xt  values  less  than  a. 

These  equations  are  difficult  to  use  for  at  least  two  reasons.  The  first  Is  that  n,  Is 
unknown  and  must  be  estimated  Indirectly  via  a.  The  second  reason  Is  Inherent  In  all 
families  of  distributions  for  which  the  density  function  can  be  Infinite  at  the  end  points 
of  the  distribution.  With  no  additional  constraints,  the  equations  always  lead  to  a  U- 
shaped  distribution  (6  >  1,  e  >  1)  with  end  points  at  the  minimum,  x^),  and  the 
maximum,  X(n^,  since  the  likelihood  function  Is  then  Infinite.  When  the  distribution  Is 
known  not  to  have  Infinite  density  at  the  lower  and/or  upper  bounds,  the  estimation 
procedure  can  be  Improved  by  specifying  the  lower  bound,  a,  and/or  the  upper  bound, 
b.  A  modified  maximum  likelihood  algorithm  Is  to  solve  the  five  nonlinear  equations 
(5),  (0),  (7),  0  —  (b-a)/(T (c/6  )1/a  +  (l-O^A)172),  and  «  =  a  +  7  fo/6  )V*  using 
Initial  estimates  for  a,  P,  7,  6,  and  c  obtained  from  one  of  the  four  sets  of  properties  In 
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Section  3.1. 


However,  the  authors’  experience  has  been  that  solving  the  system  of  nonlinear 
equations  often  results  In  numerical  complications.  Minimizing  the  average  squared 
deviation  between  the  empirical  and  fitted  distribution  function,  the  approach  taken  by 
Ozturk  and  Dale  (1985)  for  the  generalized  lambda  distribution.  Is  an  approach  worth 
Investigating.  The  authors  have  obtained  good  fits  by  repeatedly  using  one  of  the  four 
techniques  of  Section  3.1  for  several  Input  sets  of  values  and  choosing  the  parameters 
with  the  minimum  associated  goodness-of-fit  value,  as  Illustrated  In  the  next  section 
with  grouped  data  and  the  x2  statistic. 

4.  NUMERICAL  EXAMPLES 

Two  numerical  examples  Illustrate  the  family  of  distributions  In  applied  settings. 

4.1  Modeling  a  Logic  Network 

In  modeling  a  logic  network,  the  propagation  delay  of  each  component  Is  of 
Interest.  The  Texas  Instruments  TTL  Data  Book  (1973)  gives  the  minimum  time, 
typical  time,  and  maximum  time  for  a  posltlve-NAND  gate  Is  2,  7,  and  15 
nannoseconds,  respectively.  In  a  particular  situation  suppose  the  model  calls  for  1%  of 
the  gates  to  exceed  the  stated  maximum  of  15  ns  and  for  the  upper  bound  to  be  20  ns. 
The  four  required  properties  (lower  bound,  mode,  upper  bound  and  (p2,  x2)),  which  do 
not  match  any  of  the  four  cases  In  Section  3.1,  are 

a  —  mode  =  7  ns 

a  -  0i(e/6)1/2  =  lowerbound  =  2  ns 
a  +  /?  (l-"y)  (5/f)l/2  =  upperbound  ==  20  ns 
and,  substituting  (p2,  x2)  =  (.99, 15ns)  Into  eq.  (2), 


15  =  a  +  0(l-^)(6/c)l/2[l  -  ((l-.09)/(l-r))e] 

Since  these  are  only  four  properties  and  five  parameters,  one  possibility  Is  to  set  6  —  e. 


The  solution  to  these  four  equations  Is  then  closed  form,  yielding 
a  =  7  ns 

/?  =  b-  a=l8ns 
n  =  0.278 

and 

6  =  e  =  ln[l  -  (x2-a)/(/?(  1-7))]  /  ln[(i-p2)/(i-or)]  =  0.2233 


4.2  Modeling  Coefficient  of  Friction  Data 

Consider  fitting  to  the  250  observed  coefficients  of  friction  In  Hahn  and  Shapiro 
(1907,  p.  219).  The  given  frequency  classes  are 


Table  1.  The  Hahn  and  Shapiro  (1867)  coefficient  of  friction  data. 


Coefficient 

of  Friction 

Percent  of 

Observations 

Less  than  0.0150 

0.4 

0.015  - 

0.020 

3.6 

0.020  - 

0.025 

12.0 

0.025  - 

0.030 

17.0 

0.030  - 

0.035 

23.2 

0.035  - 

0.040 

18.0 

0.040  - 

0.045 

11.0 

0.045  - 

0.050 

6.8 

0.050  - 

0.055 

3.0 

0.055  - 

0.000 

1.0 

0.000  or  more 

1.6 

Set  I  of  Section  3.1  was  used  based  on  forty-five  sets  of  values  for  the  all 
combinations  of  the  properties 
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(m,  pm):  (.030,  .330),  (.0325,  .452),  (.035,  .508) 

(Xj,  pt):  (.02,  .04),  (.025,  .10),  (.030,  .330) 

(x2,  p2):  (.035,  .508),  (.040,  .748),  (.045,  .804),  (.050,  .932),  (.055,  .008) 

The  smallest  x2  value,  4.387,  was  obtained  when  a  ■=  0.0,  (m,  pm)  =  (.030,  .330), 
(Xj,  pj  =  (.02,  .04),  and  (x2,  p2)  =  (.045,  .804).  The  corresponding  parameter  values 
are  a  —  .03,  /?  =  .08115,  7  =  .330,  8  =  .1905,  and  6  =  .2300.  The  corresponding 
upper  bound  Is  b  =  .07897.  As  can  be  seen  from  the  plot  of  the  empirical  distribution 
function  and  the  fitted  distribution  function  In  Figure  3,  the  fit  using  this  crude, 
nonoptlmal  method  Is  good.  The  Johnson  system  Sy  distribution  fitted  by  matching 
moments  yields  a  curve  In  Figure  3  that  Is  difficult  to  distinguish  visually  from  the 
five-parameter  curve  shown. 
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APPENDIX:  Derivation  of  Residual  Moments 

For  k  =  l,  2,  •  •  •  ,  the  kth  moment  conditional  on  {X  >  c}  Is  derived  here.  For 
any  c  G  (-  oo,  oo),  by  definition 


E{Xk  |  X  >  c,  a  =  0}  = 


/  xK 

0 


f(x)  dx  /  [1  -  pc] 


If  Pc  <  1 
If  Pc  —  1 


Substituting  for  f(x),  we  slmpllfly  the  Integral  as  follows. 

-r)  («/«)>/*  k„  .  . 

1  xK  f(x)  dx 

■'C 

=  {-  finie/6  )1/2[1  -  (P/7)5]}*  dp 

+j'{0{l-l)(6/e)l'*ll  -  ((l-p)/(l-7)n}k  dp 

Ptt 


EH)'  ( * )  (p/7)‘^[|i 


i=0 


dp 


+  (i-7)k(6/e)k/2  Jpl 


EH)1  (f)((i-p)/(i-'7)),< 

i=o 


dp 


k  (-D‘  (f) 

(--7)k  (c/5  )k/2  £  — 


i=0 


(16+1)  t** 


+  d-7)k  (6/€)k/2 


(-1)'  (*)  (l-pu)1<+1) 

S  (le+l)  (l-7)ie 


where  pc  =  F(c;  a  =  0),  pt  =  mln{<7,  pc},  and  pu  =  max{7,  pc},  as  defined  In  the  text. 
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Figure  3.  Five-parameter  curve  fitted  to  coefficient  of  friction  data  from  Hahn  and 
Shapiro  (1967). 
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